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In the past decades, the 7-decay of giant resonances has been studied using phenomenological 
models. In keeping with possible future studies performed with exotic beams, microscopically-based 
frameworks should be envisaged. In the present paper, we introduce a model which is entirely based 
on Skyrme effective interactions, and treats the ground-state decay within the fully self-consistent 
Random Phase Approximation (RPA) and the decay to low-lying states at the lowest order beyond 
RPA. The model is applied to ^"**Pb and ^°Zr, and the results are compared with experimental data. 
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I. INTRODUCTION 

Giant resonances (GRs) have been known for several 
decades to be the clear manifestation of the existence of 
nuclear collective motion. They carry definite quantum 
numbers (spatial angular momentum L, spin S, isospin 
T) and, as a rule, they exhaust a large fraction of the 
associated energy-weighted sum rule. Accordingly, the 
macroscopic picture of a giant resonance is often thought 
to be that of a coherent motion of all nuclcons. Although 
a number of experimental data and theoretical studies 
have been cumulated, as reviewed in monographic vol- 
umes P, Q , the question still exists whether we can ac- 
cess only the inclusive properties of the GRs (energy and 
fraction of energy- weighted sum rule), or more exclusive 
properties associated with the wave function of the GR. 
Ultimately, we can say that we miss an unambiguous 
confirmation of the macroscopic picture of this collective 
motion. 

Giant resonances have a finite lifetime. Being excited 
by one-body external fields, they are as a first approxima- 
tion described by coherent superpositions of 1 particle- 1 
hole (Ip-lh) configurations. The most probable damp- 
ing mechanism is their coupling to progressively more 
complicated states of 2p-2h . . . np-nh character (up to 
the eventual compound nucleus state). The associated 
contribution to the total width, the so-called spreading 
width P^, is the dominant one. The decay width associ- 
ated with the emission of one nucleon in the continuum 
(escape width, P^) is of some relevance in light nuclei 
but much less important in heavy nuclei. The 7-decay 
width is a small fraction (w 10~^) of the total width. 
Despite this, the study of the 7-decay of GRs has been 
considered a valuable tool since about 30 years 0, 01 • 

In these works, the fact that 7-decay can be a sensi- 
tive probe of the excited multipolarity, and that 7-ejectile 
coincidence measurements can improve the extraction of 
the properties of GRs, has been thoroughly discussed. 
Generally speaking, the study of the GR decay products 
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(whether particles or photons) is probably the only way 
to shed light on the microscopic properties of the states. 
To provide an example different from the standard elec- 
tric GRs discussed in Refs. 0, 3, we can add that in 
stable nuclei or in neutron-rich unstable nuclei some in- 
formation exist on the so-called low-lying or "pygmy" 
dipole states. Their nature (collective or non-collective, 
isoscalar or isovector, compressional or toroidal) is under 
strong debate. For these, as for other states, exclusive 
decay measurements would be of paramount importance 
as they could validate some theoretical picture. 

In this spirit we present here a consistent study of 
the 7-decay of giant resonances, both to the ground and 
low-lying excited states, not considering the compound 
7-decay In the past, the theoretical study of the 
7-decay of GRs has been undertaken using frameworks 
like the Nuclear Field Theory (NFT) Q or the Theory of 
Finite Fermi Systems Q- These studies have elucidated 
the basic physical mechanisms which explain the small 
7-decay probabilities and have provided results in quite 
reasonable agreement with experiment. As we discuss be- 
low, in Ref. [6| the quenching mechanisms for the decay 
of the isoscalar Giant Quadrupole Resonance (ISGQR) 
to a low-lying isoscalar states, are clearly pointed out. 
However, these studies are based on phenomenological 
models. 

After several decades, self-consistent mean-field 
(SCMF) or density functional theory (DPT) based mod- 
els have been developed, and have reached considerable 
success for the overall description of many nuclear prop- 
erties. Among these models, we can single out those 
based, respectively, on the nonrelativistic Skyrme and 
Gogny effective interactions or on covariant (or relativis- 
tic) effective Lagrangians. Years ago, some of us de- 
veloped a microscopic description of the particle decay 
of GRs based on the use of Skyrme forces [1, @ . It is 
timely to dispose of a fully microscopic description of the 
7-decay with Skyrme effective interactions and to assess 
how large predictive power it can have, and which lim- 
itations show up. A further motivation is provided by 
the recent measurements carried out at the Laboratori 
Nazionali di Legnaro (LNL) [13. 

In this work, we develop such model and we apply it to 
the 7-decay of the ISGQR in ^ospb and 9°Zr, both to the 
ground-state and to the low-lying 3~ state. The outline 
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of the paper is the following. In Section [IT] we present the 
formalism that we employ. Section |IlT] is devoted to the 
results we obtain, also comparing them with available ex- 
perimental data and with other theoretical calculations 
found in the literature. In Section IIVI we draw our con- 
clusions and eventually, in the [Appendix] we briefly give a 
guideline for the calculation of the perturbative diagrams 
needed for the decay of RPA excited states into low-lying 
collective states, as explained in Section HIl 



II. FORMALISM 




FIG. 1. Diagrams representing the decay of the vibrational 
state \nj) state to the ground-state. 



In this Section we discuss our theoretical framework. 
The transition amplitude for the emission of a photon 
of given multipolarity from the nucleus, is proportional 
to the matrix element of the electric multipole operator 
Q\^. In the long wavelength limit which is appropriate 
in our case, this latter operator takes the form 
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In this equation, the expression for the effective charge 
due to the recoil of the center of mass of the nucleus has 



been introduced (see, e.g., Ref. [11 



The gamma decay width, summed over the magnetic 
substates of the photon and of the final nuclear state, is 
then given by 



index n. Consequently, we can write 
|nJi\/) = Oi (JM) \RPA) , 
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where 0\{JM) is the creation operator for the state at 
hand, \RPA) is the RPA ground-state, A and A^ are the 
usual creation and annihilation operator of a particle- hole 
(p-h) pair coupled to JM, X and Y are the forward and 
backward RPA amplitudes, and the symbol denotes 
the time-reversal operation (see, e.g., Ref. [l^). 

At the RPA level, in the case of the decay of the state 
\nJ) to the ground-state, we obtain for the reduced ma- 
trix of Eq. ([3]) with A equal to J, 
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where the effective charge is defined in Eq. ([T|). It is 
possible to give a diagrammatic representation of the 
ground-state decay (see Fig. [T|). 
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where E is the energy of the transition and the reduced 
transition probability B associated with the above oper- 
ator Qx^ is 
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A. The 7-decay to the ground-state 

We consider in this Subsection the decay of an excited 
RPA state (that can be, e.g., a giant resonance) to the 
ground-state. We limit ourselves to spherical systems, 
and the RPA states have quantum numbers JM (we con- 
sider natural parity, or non spin-flip, states for which the 
orbital angular momentum L is the same as the total an- 
gular momentum J); in addition, they are labelled by an 



B. The 7-decay to low- lying states 

While RPA can be considered an appropriate theory 
to calculate the ground-state decay of a vibrational state, 
the same statement does not hold in the case of a decay 
between two vibrational states. The reason is that by 
construction RPA is an appropriate theory to describe 
transition amplitudes between states that differ only by 
one vibrational (phonon) state. For other processes, like 
the one at hand, the extension to a treatment beyond 
RPA is mandatory. A consistent framework which is 
available is the one provided by the Nuclear Field Theory 
(NET) [ll,[il, since this framework takes into account 
the interweaving between phonons and single-particle de- 
grees of freedom (or particle- vibration coupling, PVC), 
considered as the relevant independent building blocks of 
the low-lying spectrum of finite nuclei. In this work, we 
consider all the lowest-order contributions to the 7-decay 
between two different phonons. This amounts to writ- 
ing and evaluating all lowest-order perturbative diagrams 
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involving single-particle states and phonon states, that 
can lead from the initial to the final state by the ac- 
tion of the external electromagnetic field. The different 
degrees of freedom are coupled by particle-vibration ver- 
tices. The Nuclear Field Theory, as mentioned in the 
introduction, has been already applied to the study of 
7-decay in Ref. @. However, the main novelty of the 
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In these equations eph is equal to the difference of the 
Hartree-Fock (HF) single-particle energies ep ~ e^, and 
V is the residual particle-hole interaction: this latter is 
discussed below, together with the expression of its re- 
duced matrix elements. In all the energy denominators 
we include finite imaginary parts 77 to take into account 
the coupling to more complicated configurations not in- 
cluded in the model space. 



present work lies in the consistent use of the microscopic 
Skyrme interaction. 

The perturbative diagrams associated with the A-pole 
decay of the intitial RPA state \nJ) (at energy Ej) to the 
final state \n'J') (at energy Ejr) are shown in Fig.[2l and 
the way to evaluate them is sketched in the lAppcndlx] 
The resulting analytic expressions read 
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In all the above equations, the matrix elements of 
the operator Qx include the contribution from the nu- 
clear polarization (consequently they carry the label pol). 
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They read 



Q 



Xpol 



E 



1 



V2A + 1 



(OHgA||n'A)(^,r^'AHF||J^) 
{Ej-Ej,)-E^,+ir] 



{i\\V\\:i,n'\){n'\\\Qm 
{Ej - Ej,) + E„, + IT] 



(6) 



where |n'A) are the RPA states having multipolarity A 
(and lying at energy £"„'), while the bare operator Q\ has 
been defined in Eq. ([1]). The polarization contribution, 
that is. the second and third term in the latter equation, 



has the effect of screening partially the external field. 
In a diagrammatic way, the bare and the polarization 
contributions to Eq. ^ are drawn in Fig. [31 

It should be noted that the diagrams of Fig. [5] are re- 
lated two by two by particle-hole conjugation, so that 
[ ^A)| is the opposite of [ ^d) [ after the substitutions h' — > p' 
and h ^ p, and the same holds for the pairs [ ^fB)]-[^c) [ 
111] - Ifi)] [I?)] - [111 [p - [JT)] and [f^P^ 

As mentioned above, in the present implementation 
of the formalism we use consistently different zero-range 
Skyrme interactions. The single-particle energies e^, and 
the corresponding wavefunctions, come from the solution 
of the HF equations. The energies (and X and Y ampli- 
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tudcs) of the vibrational states are obtained through fully 
self-consistent RPA [l^. These quantities enter the re- 
duced matrix elements associated with the PVC vertices. 
The basic one, that couples the single-particle state i to 
the particle- vibration pair j plus n J, is 

(i|!T/||j>j> = V2J + TY,x;i;vjiihjp) 

ph 

+ {-y--^^+JY;^Vj{ipjh). (7) 
Vj is the particle-hole coupled matrix element 

Vj(ihjp) = 5](-)^'^'"^+^'^-™''(.7.m,j, - mj\JM) 

{m} 

X {jpmpjh - mh\JM)vihjp, (8) 

while Vihjp stands for {jirnijh'rnh\V\j jrrijjpmp) . For the 
detailed derivation of the the reduced matrix element of 
Eq. ([7]) we refer to the Appendix of Ref. [l^. In the 
[Appendlxl of the present paper we discuss the relation- 
ships between the reduced matrix element of Eq. (O and 
the other matrix elements that enter the previous formu- 
las. In our implementation, V used at the PVC vertex 
includes the tQ,ti,t2 and ^3 terms of the Skyrme force. 



III. RESULTS 

In this Section, the results obtained from our numerical 
calculations in ^°^Pb and ^°Zr are discussed. In partic- 
ular, we focus on the 7-decay width associated with 
the decay of the isoscalar giant quadrupole resonance (IS- 
GQR) either to the ground-state or to the 3~ low-lying 
state in the two systems. We have employed four differ- 
ent Slm-me forces: SLy5 [l^, SGII [3, SkP ^ and 

LNS lig. 

In all cases, we start by solving the HF equations in a 
radial mesh that extends up to 20 fm (for ^"*Pb) or 18 
fm (for ^°Zr), with a radial step of 0.1 fm. Once the HF 
solution is found, the RPA equations are solved in the 
usual matrix formulation. Vibrations (or phonons) with 
multipolarity L ranging from 1 to 3, and with natural 
parity, are calculated. The RPA model space consists of 
all the occupied states, and all the unoccupied states ly- 
ing below a cutoff energy Ec equal to 50 MeV and 40 
MeV for 208p]_, ^°Zr, respectively. The states at pos- 
itive energy are obtained by setting the system in a box, 
that is, the continuum is discrctizcd. These states have 
increasing values of the radial quantum number n, and 
are calculated for those values of I and j that are allowed 
by selection rules. With this choice of the model space 
the energy-weighted sum rules (EWSRs) satisfy the dou- 
ble commutator values at the level of about 99%; more- 
over, the energy and the fraction of EWSR of the states 
which are relevant for the following discussion are well 
converged. 



TABLE I. Energy E of the ISGQR and 7-decay width asso- 
ciated with its transition to the ground-state. The first four 
rows correspond to the present RPA calculations performed 
with different Skyrme parameter sets, for the two nuclei at 
hand. In this case, for ^"^Pb we show both the bare from 
Eq. ([2]) as well as the renormalized value which is discussed in 
the main text. The next three rows report the results of pre- 
vious theoretical calculations for ^''^Pb. In the last row 
the experimental value for ^°*Pb from Ref. 0] is displayed, 
corresponding to the direct decay. 
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[eV] 


SLy5 


12.28 


231.54 


127.58 


15.33 


211.77 


SGII 


11.72 
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182.03 


SkP 


10.28 


119.18 
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13.09 
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12.10 
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Ref. [5| 
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Ref. [6| 


11.20 




142 






Ref. [3 


10.60 




112 






Ref. [3] 


10.60 


130±40 







TABLE II. Detailed properties of the ISGQR and of the low- 
lying 3~ state in ^''^Pb. The label Exp. indicates the cor- 
responding experimental values (the italic number after the 
value is the experimental error), these values are from Ref. 0] 
for the ISGQR and from Ref. [2l[ for the 3" state. 

2+ 3" 
E [MeV] EWSR [%] E [MeV] B(E3t) [lO^c^fm*^] 



Exp. 


10.6 


90 20 


2.6145 3 


6.11 9 


SLy5 


12.28 


69.27 


3.62 


6.54 


SGII 


11.72 


72.31 


3.14 


6.58 


SkP 


10.28 


81.79 


3.29 


5.11 


LNS 


12.10 


66.98 


3.19 


5.67 



A. Ground-state decay 

We group in Table U the results obtained for the decay 
of the ISGQR to the ground-state. In general, our calcu- 
lations reproduce the experiment quite well, without any 
parameter adjustment. They tend at the same time to 
overestimate the decay width, and this is true in particu- 
lar for SLy5; however, even in this worst case, our result 
lies within 2.5tT from the experimental value. 

More importantly, this discrepancy is entirely due to 
the fact that the properties of the giant resonance (en- 
ergy and fraction of EWSR) do not fit accurately the ex- 
perimental findings. In particular, the critical quantity 
turns out to be the resonance energy, since in Eq. ^ 
the energy of the transition is raised to the fifth power: 
consequently, an increase of the energy by 1 MeV pro- 
duces an increase of the gamma decay width by about 
50% (at 10 MeV). To substantiate this point, in the last 
column of Table U we report the values obtained for the 
decay width after having rescaled the ISGQR energy to 
the experimental value (shown in Table HI)) . In particu- 
lar, in Table HIl the fraction of EWSR exhausted by the 
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ISGQR is shown. The experimental value reported here 
is from Ref. 0] and is obtained from the measurement 
of the 7-decay to the ground state. On the other hand, 
to give an idea of the experimental uncertainty on this 
observable, we can say that in Ref. Q, from direct mea- 
surements, a value ranging from 78% to 98% is found, 
depending on the background subtraction, and moreover 
in the literature several results in the range 70% — 170% 
can be found (see e.g. Ref. [1] or [2l|). This is an indi- 
cation of systematic uncertainities that include those on 
the optical potentials used in the experimental analysis. 

We can conclude that, since for all the interactions the 
experimental value of the ground-state decay width can 
be obtained simply by scaling the energy to the exper- 
imental value, it means that this kind of measurement 
is not particularly able to discriminate between models 
more than usual integral properties. 

For completeness, in Table lithe previous theoretical 
values found in literature [11-3] are listed as well. In 
Ref. , the surface coupling model (cf. Ref. ) is used 
in order to evaluate the reduced transition probability 
and the decay width. In Ref. 0, the Theory of Finite 
Fermi Systems (cf. Ref. [1^) is implemented with a sep- 
arable interaction to obtain the decay width. In Ref. [5| , 
finally, the value is estimated from the empirical energies 
and fraction of EWSR. 



B. The quadrupole strength function 

Before we apply our beyond RPA model to a detailed 
and exclusive observable such as the decay from the IS- 
GQR to the 3~ state, it is important to test that, at the 
same level of approximation, one can reproduce more 



general quantities like the strength function of the IS- 
GQR. It has been known for several decades that cou- 
pling with low-lying vibrations is the main source of the 
giant resonance width 2J]. In Ref. |22|, calculations of 



the giant resonance strength function that take into ac- 
count this coupling have been performed, based on the 
use of a phenomenological separable force in the surface 
coupling model. We perform a similar calculation here 
by using consistently the Skyrme force SLy5, as discussed 
above. 

The probability of finding the ISGQR state per unit 
energy can be written as 
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where Ai?GQR is the real part of the sum of the eight 
contributions in Eq. (|10p . while Fgqr is the imaginary 
part of the same sum. The parameter ry corresponds to 
the energy interval over which averages are taken and 
represents, in an approximate way, the coupling of the 
intermediate states to more complicated configurations. 
In our calculation we set this parameter at 1 McV. The 



diagrams 4(a) - (d) correspond to the self-energy of the 



particle (or the hole) , i.e., the processes in which the par- 
ticle or the hole reabsorbs the intermediate excitation A, 



while diagrams 4(e) - (h) arc vertex corrections which 



describe the process in which the phonon is exchanged 
between particle and hole. If J or A is a density oscilla- 
tion, the latter contributions have opposite sign with re- 
spect to the former one, regardless of the spin and isospin 
character of A or J, respectively [24| . 



The eight diagrams 4(a) 
following expressions: 
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FIG. 4. NFT diagrams contributing to the strengtli function of the giant resonance. 
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The result for the probability of finding the ISGQR, 
calculated by including in the diagrams an increasing 
number of intermediate phonons is displayed in Fig. [5] 
The RPA model space is the same used for the compu- 
tation of the decay width. Phonons with multipolarity 
ranging from to 4 and with natural parity (—1)^, have 
been considered. Only those having energy smaller than 
30 MeV and fraction of the total isoscalar and isovector 
EWSR larger than 5% have been selected as interme- 
diate states. The most important contribution to the 



spreading width F"^ of the resonance is given by the low- 
lying 3~ state, while the other phonons do contribute 
basically only to the energy shift. We obtain eventually 
a spreading width F^ of the order of 2 MeV, and the 
energy centroid of the resonance is shifted down, as com- 
pared to the RPA value, to 10.9 MeV. These results are 
in good agreement with the experimental findings, that 
give a spreading width of 2.4 ± 0.4 p5j |. 
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TABLE IV. The various quenching factors that combine to 
produce the decay width from a typical particle-hole dipole 
transition, for ^"^Pb. The decay width reported here refers to 
a cutoff of 5% on the percentage of isovector EWSR for the 
dipole states. The same quantities from @] are displayed. 





SLy5 


SGll 


SkP 


LNS 


Ref. [6| 


Ph transition [eV] 


10=* 


10=^ 


10=* 


10=* 


10=^ 


Recoupling coefficient 


3 


3 


3 


3 


3 


■K - V cancellation 


5 


4 


3-4 


4 


4 


p - h cancellation 


3-4 


2-3 


2-3 


3-4 


2-3 


Polarization 


6 


3 


7-8 


4 


15 


[eV] 


3.39 


29.18 


8.34 


39.87 


3.50 



FIG. 5. (color online) Probability P per unit energy to find 
the ISGQR at an energy E in '^"^Pb. Each line corresponds to 
the probability obtained when the phonons listed in the legend 
are used as intermediate states (for the selection criteria of 
the phonons we refer to the main text). The red-solid line is 
the probability we get when only the low-lying 3^ phonon at 
3.61 MeV is considered as intermediate state. The label RPA 
(black-dashed line) refers to the RPA result, in which none of 
the diagrams in Fig. |4]are taken into account, but a lorentzian 
averaging with functions having 1 MeV width is introduced. 



TABLE III. Decay width to the low-lying 3 for the inter- 
actions used, calculated including beyond RPA contributions 
for the two nuclei ^°^Pb and 9°Zr. In particular, for ^°*Pb 
the results from Ref. 0] and Ref. are also listed and in the 
last row, and the experimental value from Ref. [^1 is provided 
as well. 





208pb 




90 Zr 




Interaction 


-Etrans [MeV] 


r, [eV] 


Strans [MeV] 


[eV] 


SLy5 


8.66 


3.39 


12.51 


5.81 


SGII 


8.58 


29.18 


12.16 


50.58 


SkP 


6.99 


8.34 


10.42 


5.14 


LNS 


8.90 


39.87 


12.72 


16.95 


Ref. [6| 


8.59 


3.5 






Ref. [7| 


7.99 


4 






Ref. |3| 


7.99 


5±5 







C. ISGQR decay to the low-lying 3 state 

In Table IIIII the results obtained for the decay of the 
ISGQR to the low-lying octupole state in ^''^Pb and ^°Zr 
are shown. These correspond to a choice of the lower cut- 
off of 5% on the dipole EWSR of the states considered 
to calculate the polarization and a parameter rj = 2 MeV 
(for 208pb) and 77 = 3.5 MeV for ^"Zr; these inputs will be 
clarified and discussed later in the text. For complete- 
ness, the values found in literature @, 0| for ^°®Pb are 
listed as well. These are all theoretical results obtained 
using different models: in Ref. 0], the Theory of Finite 
Fermi Systems (cf. Ref. [2^) with a phenomenological 



interaction is used to calculate the decay width, while in 
Ref. the decay width is obtained by means of the NFT, 
with a separable interaction at the particle- vibration ver- 
tex. We discuss here in detail the results we obtain for 
the decay of the ISGQR in ^°®Pb. It can be noticed that 
only two interactions, namely SLy5 and SkP, can rea- 
sonably reproduce the experimental value for the decay 
width. Nevertheless, all these forces are able to produce 
a total r.y (ISGQR — !• 3^) which is only a few percent of 

(ISGQR — > g.s.), as the experiment indicates. 

In order to understand which of the factors that ap- 
pear in the several contributions to the decay width has 
a major effect on the resulting values, we have analyzed 
the sensitivity to the physical inputs in great detail. Ta- 
ble |IV] displays, for the four forces used, the contribution 
of the several factors included in Eq. ([5]). Similar fac- 
tors from Ref. @ are provided as well. The decay width 
that is obtained considering a typical particle-hole tran- 
sition is of the order of sa keV and can be qualitatively 
accounted by means of the Weisskopf estimation for the 
reduced transition probability of a single particle excita- 
tion [27} . The label recoupling coefficient indicates the 
quenching deriving from the mismatch of angular mo- 
menta of the particles involved in the process. Then, 
because of the isovector nature of the operator ([1]), dia- 
grams involving protons and neutrons have opposite sign 
and partially cancel each other. Moreover, diagrams in 
which the operator acts on a particle line must have an 
opposite sign to the ones in which it acts on a hole line, 
reflecting the correlations between particles and holes in 
vibrations [24j . resulting in a compensation of the two 
contributions. Eventually, the polarization contribution 
deriving from the screening of the external field by 
the mediation of the giant dipole resonance, represent 
a further and more important quenching of the original 
decay width, giving then a final width of the order of 
electronvolts. 

We have studied in particular which assumptions and 
choices affect the quenching associated with the polariza- 
tion contribution. First of all, in Table |Vl the variation 
of the 7-decay width r..^ with the parameter rj = 



9 



TABLE V. The effect on the 7-decay width T.y of the width Td of the intermetiate dipole states. AU the states exhausting 
the EWSR for more than 5% are considered. The decay width is an almost monotonically non-decreasing function of this 
parameter, as expected from the Bohr-Mottelson model [261 ]. 





Td [MeV] 


0.01 


0.1 


0.5 


1.0 


2.0 


3.0 


4.0 


5.0 




SLy5 


2.34 


2.35 


2.46 


2.69 


3.39 


4.32 


5.34 


6.36 




SGII 


27.28 


27.19 


27.00 


27.25 


29.18 


32.68 


37.21 


42.25 


r, [eV] 


SkP 


7.35 


7.35 


7.40 


7.59 


8.34 


9.54 


11.08 


12.86 




LNS 


37.93 


37.82 


37.55 


37.76 


39.87 


43.84 


49.03 


54.86 




IMeV] 



is known to be not reliable for non-collective states, and 
secondly, introducing them would oblige to take into ac- 
count the issue of the Pauli principle correction. We 
then choose 5% as lower bound of the isovector and 
isoscalar EWSR, in keeping with several previous works, 
e.g., Ref. [13. 

A similar analysis carried out on the 7-dccay of the 
ISGQR in ^°Zr into the lowest 3~ state would bring to 
analogous conclusions: the most important effect is the 
polarization of the nuclear medium through the excita- 
tion of dipole states. Even in this case the general result 
is that the decay width to the octupole state is few per- 
cent of the one to the ground-state. Results from the 
previously mentioned recent experiment [lol | , are not yet 
available. 



FIG. 6. (color online) Polarization contribution in the separa- 
ble framework as a function of the parameter Td (solid line). 
The arrow indicates the value used in Ref. 0| . The points are 
the analogous factors obtained within our model in which all 
the dipole states having a fraction of EWSR larger than 5% 
are considered and each one is given a width Fd = 2 MeV. 



that appears in Eq. © as imaginary part of the energy 
denominator, is discussed. If only a single dipole inter- 
mediate state is considered, as in Ref. |6|, this parame- 
ter should be set equal to the IVGDR width (-4 MeV); 
since in our model, the dipole strength is fragmented, we 
should take a smaller value and we give here the trend 
of the decay width as a function of this parameter. As 
indicated by the plot in Fig. [51 the polarization factor 
(and consequently the decay width) should be monoton- 
ically non-decreasing when Tjj increases and reaches a 
roughly constant value as F^) goes to zero. In the same 
plot, the points represent the polarization factors that 
we obtain using the value 2 MeV for the parameter F^, 
but including all the dipole states having a fraction of 
EWSR larger than 5%. This value has been chosen in 
order give a width of the RPA dipole states, each convo- 
luted with a lorentzian of width equal to F^j, similar to 
the experimental IVGDR width. The polarization that 
we get is then consistent with the one of Bohr-Mottelson 
model [2^ , indicated with the arrow in Fig. [51 

We need a lower cutoff on the collectivity of the in- 
termediate states for at least two reasons: firstly, RPA 



IV. CONCLUSIONS 

Our work is motivated by the fact that we deem it 
is timely to dispose of a fully-microscopic description of 
some exclusive properties of giant resonances, like the 
7-decay. In particular, the 7-decay has been studied 
in the past decades using only phenomenological mod- 
els. Therefore, we have implemented a scheme in which 
the single particle states are obtained within HE, the vi- 
brations are calculated using fully self-consistent RPA 
and the whole Skyrme force is employed at the particle- 
vibration vertices. We treat the ground-state decay 
within the fully self-consistent RPA and the decay to 
low-lying collective vibrations at the lowest contributing 
order of perturbation theory beyond RPA. 

We have applied our model to the 7-decay of the 
isoscalar giant quadrupole resonance in ^°^Pb and ^°Zr 
into the ground-state and the first low-lying octupole vi- 
bration. In particular, in ^"^Pb, in the case of the ground- 
state decay, we find that our outcomes are consistent with 
previous theoretical calculations, based on phenomeno- 
logical models, and with the experimental data. In par- 
ticular, all the Skyrme parametrizations give a 7-decay 
width to the ground-state of the order of hundreds of elec- 
tronvolts, though, at the same time, they tend to overes- 
timate it: these discrepancies are due to the fact that the 
energy of the resonance does not completely agree with 
the experimental data. For this reason, we conclude that 
the 7-decay to the ground-state is not so able to discrim- 
inate between different models, at least not more than 
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other inclusive observable (as energy and strength). 

On the other hand, the 7-decay to low-lying collective 
states is more sensitive to the interaction used. As a mat- 
ter of fact, only two interactions (namely SLy5 and SkP) 
manage to achieve a decay width of few electronvolts, 
consistently with the experimental finding. In the case of 
SLy5, this fact is consistent with the good features that 
this parameter set has, as far as spin-independent pro- 
cesses are concerned (correct value of the nuclear incom- 
pressibility, reasonable fit of the neutron matter equa- 
tion of state, good isovector properties). Nonetheless, the 
other interactions give a width that is of the order of 
tens of electronvolts and it is very much quenched with 
respect to the decay width associated with a single parti- 
cle transition. It is quite remarkable that our calculation, 
being parameter-free, reproduces numbers that are sev- 
eral orders of magnitude smaller than the nuclear scale of 
« MeV. In particular, the description of the dipole spec- 
trum is a crucial point, because small differences in the 
strength of the dipole states, introduced as intermediate 
states, change significantly the polarization of the nuclear 
medium. For ^^Zr, the general conclusion is similar: the 
7-decay to low-lying collective states seems to be a good 
observable to test the quality of different Skyrme models, 
being very sensitive to the description of the polarization 
of the nuclear medium. 



Appendix: Calculation of the diagrams associated 
with the decay between vibrational states 

In this Appendix, we provide some details about the 
calculation of the diagrams shown in Fig. [21 

Within the PVC theory, four particlc-phonon ver- 
tices are possible (cf. Fig. [7]), depending whether the 
fermionic states involved are particles or holes. They 
are related by the particle-hole coniugation operator (cf. 
Ref. [2^ ) , so that all the vertices can be brought back to 
Ap, From the Appendix of Ref. [ll], we get for the first 



vertex A^,p 

Ap'p = {jp'mp>\V\jpmp,nJM) 

(A.l) 

For example, let us now consider the vertex A^^, . We 
can move the hole states from the initial to final state 
(and viceversa) by adding an appropriate phase factor. 

Kw = {Uh'm.h')-^\V\{jHm,,)-\nJM) 

(A.2) 

which is equal to A^/^ after the identification h' ^ p and 
h ^ p' , except for the phase factor. Similar relations can 
be established between Ap,p and A^^ or A'f^p. Moreover, 
the vertices in which the phonon is created instead of 
annihilated can be derived from these latter by adding a 
phase factor (— )'^"'"*^, changing the sign of the projection 
M of the angular momentum J of the phonon and using 
the following expression for the reduced matrix element 
of the interaction 

{^,nJ\\V\\J) = V2jTT{-y+'^-''Y.X;i!Vj{jMp) 

ph 

+{-y'^-^.+JY;^fVj{jpzh). 

These PVC vertices are then used to evaluate the di- 
agrams in Fig. [21 In the following, one of them (namely 
Fig. [ ^e)| is calculated in detail. For each particle- 
phonon vertex, we have a reduced matrix element of the 
interaction multiplied by a 3-j symbol that takes care of 
the coupling of angular momenta. Moreover, the single- 
particle operator Q\ brings another 3-j symbol and a ma- 
trix element. Eventually, the last 3-j symbol, matching 
the angular momentum of the initial state, of the final 
one and of the operator comes from the Wigner-Eckart 
theorem, since we need a reduced matrix element. The 
energy denominators are obtained by using the rules of 
second-order perturbation theory. 
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FIG. 7. Particle-vibration coupling vertices 



Jh Jp' A\ V I ^hp' 

ruh nip' nj (Ej - eph+iT]){Ej -hujj, -ep'h + iif) 



(A.3) 



r 



The four 3-j symbol can be summed in one 6-j symbol by usual relations (see e.g. Ref. j28| ) 
I 



h-mh f jp' jp J' \ f jp jh J\ f jh jp' A\ / J A J 



MpM' 



-nipi Trip —A/' 



mil M J ynifi mp> /i J ii — il/' 



J A J 

v'^TM^ \jp' Jp jh 



(A.4) 



We then finally get Eq. ()5.eI 
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